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Abstract The Convolution and Master equations governing the time behavior of the 
term structure of Interest Rates are particularly simple for continuous variables. The 
infinitesimal generator of the generalized Markov process is usually a distribution. 
We believe that the discretised forms of the equation and of the generator are better 
suited to compare with actual distributions. They help to avoid the Gaussian-like 
tail behavior generally derived form the continuous equations with a finite number of 
Markov processes. In this paper, the notion of discretised Seed is introduced which 
naturally leads to an infinite superposition of Markov processes and hence allows a 
power (rather than exponential) decrease of the empirical probabilities of the varia- 
tions of the interest rates. The discretised theoretical distributions of probabilities 
matching the empirical data from the Federal Reserve System (FRS) are deduced 
from a discretised seed which enjoys remarkable scaling laws. In particular the tails 
of the distributions are very well reproduced. These results may be used to develop 
new methods for the computation of the value-at-risk and fixed-income derivative 
pricing and suggest the appearance of the critical exponents related to models based 
on self-organized systems. 
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1 Introduction 



The accuracy of the interest rates variations modelling is an important issue espe- 
cially in the context of the evaluations of the value-at- risk and marked-to-market 
positions in trading floors. 

In two recent articles pQ, j2], it was shown using empirical data published by 
the governors of the Federal Reserve System [3j from 1962 until 2002 that the term 
structure of interest rates decreases essentially for large variations of the interest 
rates as a power and that this power is of the order three to four. Moreover the 
distributions seem to obey simple approximate scaling laws as functions of the initial 
interest rate, of the lag and of the maturity. These findings invalidate many models 
which predict distributions having either very short tails, generally exponentially 
decreasing, or very long tails as do Levy type structures [S], [I], |E] 

In this paper, a theoretical model is built to serve as a basis for computing 
the distribution of the variation of interest rates in terms of a few fundamental 
"microscopic parameters" whose meaning will be highlighted. At the basis of the 
theory, the notion of "seed" (a discretised form of the infinitesimal generator ^]) 
is introduced. It is closely related to the variation of the interest rates for a very 
short but finite time intervals. We believe that this discretisation is an important 
issue, closely linked to the fact that empirical interest rates are, often by regulation, 
expressed as integer multiples of the basis point and obtained by averages performed 
on a finite set of discrete times. 

In order to simplify the presentation and highlight preceding results available 
in the litterature jl], the problem is set in terms of continuous variables. Later, 
to allow numerical simulations and come closer to the empirical distributions, the 
variables are discretised. 

All along this article, in order to make a connection to a real situation, we have 
chosen to refer systematically to an application of our ideas to the FRS data [3]. 
Needless to say, we expect our analysis to be extendable, mutatis mutandis, to many 
other situations. 

2 The basic equations with continuous variables 

In this section, the basic equations, which govern the continuous time propagation 
of the term structure of interest rates, are briefly reviewed, analysed and commented 
from a eco no-physicist's point of view. 

Suppose that, the interest rate for a certain maturity [m] has the value J at 
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time t - We want to study the normalized density of probability (as seen at time t) 

P [ r\t f ,i f ,t ,i ) (i) 

that, at a later (final) time £/, thus after a lag 

L = t f -t , (2) 

the interest rate has a value If. In principle, this density of probability p t is an 
unknown function of the five "continuous variables" t,tf,If,to,Io. In order to sim- 
plify the notation we will restrict ourselves to a given maturity and suppress the 
corresponding upper index. 

The compounded probability (at time t) Pt,t f ,i a <if<h,to,io that, starting with an 
initial interest rate Iq at the initial time to, the final interest rate If at the final time 
tf is in the interval between I a and 1^ is given by the integral on this interval of the 
probability density 

Pt,t f ,i a <i f <i b ,to,io= Pt(tf,If,t ,I )dIf (3) 
•''it 

with, obviously, by normalization 

Pt,tf,-oo<If<+oo,t ,Io — 1 ■ (4) 

We will now make the hypothesis that the interest rate variations satisfy some 
market laws which are rather stable and are governed by sufficiently smooth equa- 
tions. Let us try to state more precisely and justify our simplifying assumptions. 

2.1 Time translation invariance 

This is a very delicate hypothesis. It is equivalent in saying that, whatever be the 
political or economical situation the average behavior will be identical. The effects 
of the exceptional situations (political or economical crisis) which usually may lead 
to seemingly incoherent variations of the interest rates are accurately taken into 
account by the tails of the functions which are used. The postulate is that, even if 
there are extreme situations of various importance, all in all they connect smoothly 
with the more normal situations which prevail during the peaceful times. It is 
precisely these extreme situations which are the prime reason for the fat tails of the 
distributions. And they must be incorporated correctly by the model. To repeat, 
the hypothesis is that what happens in the exceptional situations is well taken into 



2 



account by the smooth tails of the distributions. There is no abrupt transition 
between really exceptional situations and what we would call the normal situations. 
Between the extremes, there are situations of intermediate seriousness which lead 
for the distribution to a smooth passage from a restful period to a chaotic one. But 
let us make our arguments more precise. 

Mathematically, for any of the times t, t , tf, the global time translation operation 
is given by 

t' = t + T (5) 
f = t + T (6) 
t' f = t f + T. (7) 

where T is the value of the translation time (minutes or days or months or years 
later). Suppose that at some time t' later than t given by the translation time T 
the interest rate I' is again exactly equal to 7 . Technically, the time translation 
invariance demands that the densities of probability to obtain the same final interest 
rate I'j = If be equal 

Pi/ (t'f, J f, t'o, I ) = Pt(t f , If, t , Iq) (8) 

whatever be T. The financial interpretation of this equality is that the market laws 
are stable in time. If the same situation is reproduced later, it will evolve with the 
same probabilities. 

It is not difficult to prove mathematically that this condition essentially implies 
that the density of probability does not depend on the three time variables t,t ,tf 
separately but only on two of their differences, say t — t and the lag 

L = t f -t . (9) 

Let us finally discuss the dependence in t which is the time at which the information 
is needed. Obviously if t is larger than to, i n principle all the information of the what 
the actual interest rates have been between to and t is available and known. The 
related probabilities are nor useful any more. If t is smaller than to it means that 
we want to evaluate the probabilities for evolutions of the interest rates happening 
at later times. Again, if the market laws are stable, the evolution during a fixed lag 
L as seen at the time(t = to) or at a earlier time (t < to) should be identical. The 
same initial state (I at to) should lead statistically to the same final state (// at 
tf). This means invariance under a further logically different time transformation 



3 



((HI) while t and tf are kept fixed 

t' = t + T (10) 

t' = t (11) 

t' f = t f (12) 

This implies that the density of probability is independent of its index t. 

Gathering the result of the two invariances, we find the time invariant restricted 
form of the density 

p t (t f , I f , t , 1 ) = p(t f - t , If, I ) . (13) 

In other words, it does not depend independently on t, tf and t but on the lag only. 

The invariance under the two continuous time translations thus implies that the 
density of probability pt is reduced to a function of three continuous variables only: 
the lag, the final interest rates // and initial interest rate Io- The variation V of the 
interest rate during the lag L is defined as 

V = If — Iq . (14) 

For later convenience, a last variable change is performed to define the time transla- 
tion invariant density of probability p as a function of the lag, of the initial interest 
rate and of the variation of the interest rate. From now on, the form 

p(L, V, Io) =p(t f - h, V + Io, Io) (15) 

will be used as the basic interest rate distributions. 

2.2 The normalization of the probability 

The total probability (@J that after the lag L the rate If has any value must be 
equal to one. This implies the normalization 

/+oo r+oo 
p(L, If - Io, Io) dlf = / p(L, V, Io) dV = l. (16) 
-oo J — oo 

This normalization should hold whatever be the lag and whatever be the initial rate. 

2.3 The composition of the probabilities 

The basic equations, which govern the composition of the probability densities p 
(essentially the Chapman-Kolmogoroff equations 4j), are well-known. Let us write 
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them in our notation. Intuitively, consider three times to,ti,tf where the interme- 
diate time ti lies between the initial time to and the final time tf. 

t <U<t f . (17) 

The initial rate is Iq. 

During the time lag L\ = (ti — to) the interest rate has a density of probability 
p{ti — to, Ii — Jo, Iq) to reach the intermediate value Ij. Then starting from the 
intermediate time ti up to final time tf (i.e. during the second time lag L 2 = (tf—ti)) 
the intermediate observed interest rate lj at time ti has a density of probability 
p(tf — ti, If — Ii, Ii) to become If at the final time tf. Since the interest rate at the 
intermediate time ti can take any value, the density of probability starting from the 
rate Io at initial time to to end up with a rate If at the final time tf is given by the 
integration on Ii at the intermediate time (convolution of the probabilities) 

/+oo 
p{t f - U, If - Ii, Ii) p(ti - t , Ii - I , 1 ) dli . (18) 

This is the basic equation which the probability distribution has to fulfil. It shows 
how the probability distributions of two successive lags L\ = (ti — to) and L 2 = 
(tf—ti) compose to form the probability distribution for the lag L = tf—t = Li+L 2 . 
Eq. (jl8|) should hold whatever be the intermediate time. 

It is convenient to rewrite the equation by using the new variables I, V and the 
new integration variable W 

I = 
V = 
w = 

In = 

as 

/+oo 
p(L 2 ,V -W,I + W)p(L 1 ,W,I)dW (20) 
-oo 

This is the basic equation. 
2.4 Initial conditions 

The probability distribution has to satisfy an initial condition which can be described 
in terms of the Dirac 5 distribution (see Appendix A). Indeed if the initial rate at 
time to has the value Iq (whatever I is), at tf = to (i.e. after a zero lag L = tf —t = 



If 



ti — to 



If 
h 
L 2 



I 

V + I 
W + I 

tf - ti 



(19) 
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0) we know for sure that the rate is still Jo- The density of probability for the final 
rate to be different from J is zero. In other words, when the lag is zero, the rate 
has not moved. It is still Jo with probability one. 
Mathematically, this implies that 

p(0,V,I ) = 5(V) . (21) 

Some properties of the distribution S(V) are given in Appendix (|2J together with 
a few useful approximations in terms of more conventional functions which will be 
used later. In particular, the Dirac distribution is normalized 

/+oo r+oo 
p(0, V, J ) dV = 5{V) dV = 1 . (22) 
-oo J — oo 

As it should, in agreement with the composition of probabilities (fTHjl evaluated either 
for ti = to or for ti = tf, the Dirac distribution, composed with any distribution, 
satisfies the identities 

/ + O0 
p(tf - t Q , If - It, li) p(0, li - J , J ) dli 
-oo 

/+oo 
p(0, If - I h h) p(t f - to, U - J , Jo) dli (23) 
-oo 

showing again that an evolution during a zero lag is, in fact, not an evolution. 



2.5 The seed 

By using the convolution equation JTBJ), if one knows the probability distributions 
p(e, V, J ) for a given lag L = e, whatever be the value of e, the probability distri- 
bution for p(2e, V, Jo) for a lag of 2e can easily be computed by the convolution of 
p(e, V, Jo) with itself. Then p(3e, V, Jo) is obtained by the convolution of p(2e) with 
p(e, V, Jo). By successive iterations p(ne, V, Jo) can be computed for any positive 
integer n. Hence, if one knew the distribution of probability for a very small lag, e 
much smaller than the empirical lags, the distribution for a lag L much larger than 
e could be computed approximately by simple successive integration. Essentially a 
number of integration equal to the integer n closest to (L/e—1). This approximation 
would become better and better when e — > 0. 

It is then tempting to let e go to zero and to take the first order approximation, 
which we call the seed S(V, J ) for reasons to follow, as 

styM = lim p(^,w-K<w.) 

e^O e 

= d L p(L,V,I ) . (24) 

L=0 
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Since the V^-integration of the two p terms in the right hand side are equal (before 
e is put to zero) by (JTBJ) , the seed S satisfies the normalization condition 



This seed is essentially related to what is often called the "infinitesimal genera- 
tor" of a Markov process in analogy with the infinitesimal generator familiar in the 
theory of Lie groups and semi-groups. It allows, in principal, the computation of the 
time evolution for a finite time from the knowledge of what has happened during an 
infinitesimal small amount of time. 

We have however decided to call it seed for two reasons. First, all the variables 
(interest rates and times) are discretised. Thus, the meaning of the seed does not 
coincide with the meaning of the infinitesimal generator as it does not correspond 
to an infinitesimal time anymore. Second, Markov processes are usually based on 
random walk models. For a finite superposition of Markov processes (arbitrary large 
but finite) the distribution of probability turns out to be a finite superposition of 
exponentials (short tails) (see ©)• Our generalisation will in fact be equivalent to 
an infinite superposition and hence will transcend exponentials and hence escape 
the ill-fated short tails. 

If a Taylor expansion is used, in first approximation, for L = e very small, the 
probability distribution p can written 



In the continuous case, this is a distribution. Nevertheless in the discretised case, this 
form, with function rather distributions, can be used as a very good approximation. 

In the rest of the section, we discuss a few facts related to the presence of 
distributions rather than functions. The reader which is not interested in these 
mathematical considerations should jump immediately to section (J3J). Indeed, tech- 
nically, it should be stressed that, as the initial condition (12 lj) is a distribution rather 
than a function, we should expect that the seed is also a distribution with support 
restricted to V = 0. 

2.6 The Master Equation. Continuous variables 

If the seed is known as a distribution, the density of probability of the interest rate 
variation can easily be computed by integrating the master equation as an integro- 
differential equation 




(25) 




(26) 



J —oo 



)S(Ii-I ,I )dIi 



(27) 
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which follows from the convolution equation (|18|) . With the change of variables (|19|) 
the master equation can equivalently be written 



Technically, one obtains Eq. (j27|) by differentiating the convolution equation (fTHj) 
with respect to t« and letting tj — > t . The master equation (}2"T|) is a consequence of 
the convolution equation. Conversely, in Appendix JBJ), it is shown, formally, that 
the solutions of the master equation satisfy the convolution equation (|18j) . 

3 The Gauss distribution as the solution for the 
simplest seed 

The simplest possible (joTj) (see Appendix (jXJ)) seed is the distribution 



This can be justified naively as follows. During a very small amount of time e one 
expects the variation of the interest rate V to be very small, say at most of the 
order p where p decreases and goes to zero with e. Introduce the three intervals of 
length p, the left interval L = [— 3p/2, — p/2], the center one C = [—p/2, +p/2] and 
the right one R = [+p/2, +3p/2]. Outside these three intervals the probability is 
essentially zero as is the density of probability. In the left and right intervals the 
probability is small say of the order k and in the center it must be 1 — 2k to conserve 
the normalization of the probability. If the density of probability is then supposed 
to be constant in the three intervals, the seed becomes n/p in the left and right 
intervals and —2n/p in the center interval. Letting the time interval e go to zero 
induces p to go to zero and the limiting seed becomes, up to k, the second derivative 
of the 5 distribution. A Markov process is a finite superposition of terms of the form 
(j2U|) with possible higher (even) derivatives of the 5 distribution. 

Taking this simplest form of the seed, the integro-differential equation (J27j) be- 
comes a much simpler partial differential equation 




(28) 



S(V,I ) = n(I )d 2 v 5(V) . 



(29) 



d L p(L, I f - I , 1 ) = k(I ) d 2 Io p{L, I f - I , 1 ) 



(30) 



which is also written 



d L p(L, V, Jo) = «(/ ) (d 2 v p(L, V, J ) - 2 d v d Io p(L, V, I ) + d\ 



p(L,VJo)) ■ (31) 
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This equation is still rather complicated. 

It is shown in Appendix (|Cj). assuming no Jo dependence and simplified scaling 
laws (JHHj) and (|51jl . that the solution of the master equation for the seed (}2T?j) is of 

Gaussian type 

P(V) = ^J==e-£ . (32) 

This solution as highlighted in [I] and [2] is completely excluded by the obser- 
vations as are all Markov's approaches which lead to finite superpositions of such 
Gaussians. 



4 Some considerations about the discretisation of 
the problem. Numerical integration 

4.1 Why discretisation? 

The choice of discretising the distributions is not done here for convennience but 
is motivated by the discreet nature of the available data and the way they are 
produced by the laws of supply and demand. Naturally, market makers are always 
working in units of "something, 1/32 of $, basis point . . .. Thus, the natural variation 
measure units are by nature discreet. The same discussion applies to Federal Reserve 
System data In this case, we are quite naturally led to discretise both the time 
variable and the interest rates variable. As a consequence the master equation be 
better discretised Indeed, first, by regulation, the interest rates are expressed by 
integers in basis point (10~ 2 percent) and hence naturally discreet. Then, usually, 
the determination of the daily average interest rates is compounded by a finite 
set of individual contracts finalized at certain times. These times do not occur at 
fixed intervals. But they are clearly not continuous and probably better taken into 
account by a set of discrete times at regular intervals rather than by a continuous 
time. This is a critical point. We show that the discretisation leads to extremely 
accurate results. 

Moreover, technically, given a seed which is a distribution, the convolution equa- 
tion ([2 L)|) or the master equation ()28jl are usually very difficult or impossible to solve 
analytically. Instead, in this present work, a numerical procedure has been imple- 
mented. This numerical integration depends on the discretisation of the continuous 
equation using discretised variables on a grid. 

We have focussed our attention to the following specific form of the continuous 
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equations (|20|) for L = L\ = L 2 

/+oo 
p{L, V - W, I + W) p(L, W, I) dW . (33) 
-oo 

In the following, the ideas and the delicate points are briefly described. The pre- 
cise application of the formalism to the FRS data jH] and the results are summarized 
and highlighted in section ((Zj). 

In this section are presented in turn the ideas which are needed to perform the 
approximate numerical integration of the equation. The continuous equations are 
discretised by associating a finite lattice (a grid) to the continuous variables. The 
initial conditions and the seed are then defined on a finite set of points at the nodes 
of the grid. This allows then to define the x-squared merit function. 

4.2 The ideas behind the grid 

The equation (|33J) depends on three continuous variables. In order to perform nu- 
merical integrations, the three variables have to be discretised. 

The Lag grid 

Suppose that the distribution is known for some value sl of the lag L. From 
the first application of (}3*3*j) one obtains the distribution for L = 2s l- Applying 
(jSni) a second time allows the determination of the distribution for 4sl = 2 2 sl- 
The iteration of the procedure produces the distribution for any L of the form 
2 1 sl i.e. on a the L-grid of values 

L = 2 l -s L l = 0,...,N L . (34) 

Our choice of the Lag grid 

We have chosen to obtain a lag L = 1 day after Nl = 10 iteration. This leads 
to a value of ^ 

i.e. about 30 seconds. This is a perfect choice if individual contracts are 
concluded at this rate. If it turns out that the lapse between two successive 
contracts is smaller, Nl should be increased. 

The V grid 
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For the V variable, the best choice is a regularly spaced grid consisting of 
points of the form 

V = v ■ s v (36) 

where sy is the step size in V and v is an integer. A smaller sl should lead to 
a choice of a smaller sy- 

Our choice of the V grid 

We have chosen to take s v to be 

/ 1 \ 10 

sy = \ — \ basis point (37) 

i.e. around (775) 10 basis point. Indeed a larger sy turned out not to be small 
enough while a smaller one was not necessary and increased the computation 
time and memory requirements with little reward. 

The Iq grid 

For the Iq variable, the best choice is a regularly spaced grid consisting of 
points of the form 

Io = i-S! (38) 
where sj is the step size in I and i is an integer. 

Our choice of the J grid 

Since all the parameters vary with I rather weakly, a choice of 

si = one percent (39) 

is quite sufficient. 

In our actual computations due memory size problems, we have been led to adapt 
the sy step size to the level / of the iteration. The step size is progressively increased 
in such a way as to become one basis point at the tenth iteration. This required 
delicate numerical adjustments 
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4.3 The ideas behind the discretised equations 

Obviously, on the grid (J53j) . (jHSJ), (JHHj) the equation (J3~3~j) becomes 

{oo 
V (2's L , (v - w)s Vl isj + ws v ^j x 
w=— oo 

x p^s^wsv^s^} . (40) 

To be on the grid, the argument i sj + w sy which appears in the first p on the right 
hand side must be an integer multiple of s/. This would at first sight imply that sy 
is an integer multiple of s/. In fact, our numerical implementation uses, for technical 
reasons, a more elaborate grid than (JHUj) in the V and a better adapted form of the 
convolution ()40|) together with smoothing and interpolation techniques. Moreover 
the normalization implies 

oo 

]T p(2 l s L ,vs v ,i Sl ) = l. (41) 

v=— oo 

4.4 The ideas behind the discretised initial condition 

The initial condition (|2~T|) on the grid is taken as a step-type f!69|) approximation of 
the Dirac 5 distribution, as explained in Appendix Namely 

1 

p(0,0,is/) = — , 
sv 

p(0,vs v ,isi) = forv^O. (42) 

4.5 The ideas behind the discretised seed 

Using the discretisation initial condition, for e = sl chosen sufficiently small, the 
Taylor expansion (J26j) of p to the first order in Sl 

p(s L , v s v , i si) = p(0, v s v , % sj) + s L S(v s v , i s/) . (43) 

The resulting p(sl, v sy, i si) must be a true function (not a distribution anymore). 
It must be positive and normalized (}4~Tj) 

oo 

£ S(vs v ,i Sl )=0. (44) 
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Hence, S(v sy,i sj) is expected to be a step-type function of vsy with the following 
properties 

S(0, is/) = -**r,S(vs v ,isi) 

S(vs v ,is!) > forw^O (45) 

S{0,t Sl ) < — . 

slSv 

The form of the seed inferred by the FRS data is presented and discussed in sec- 
tion ©. 

4.6 The ideas behind the x-squared function 

In order to determine the seed free parameters, a x _sc l uare d minimization method 
is used. 

As usual the generic x-squared function is defined by 

x 2 = y (AWoryQ) ~ jVdata(fc)) 2 ^ 

where the sum is performed on the generic discrete label k indexing all the available 
data. In (|4b|). iV data (/c) and N theory (k) are respectively the number of events observed 
and predicted for the label k. 

The natural error (a (fc)naturai) is usually taken as if the distribution was of Poisson 
type i.e. 



0"(/c) nat ural = \J N data (k) if iVdataO) 7^ 

^(AOnatural = 1 if ^data(*0 = ■ (47) 

The choice of the error function requires more attention in the present case. As 
has been already stressed, the theory should take into account and reproduce as 
accurately as possible the tails of the experimental distributions. This is achieved 
by choosing in the x-squared a modified form of a{k) 



O"(£0modificd = (N A&ta {k)N theovy {k)) l,4: if iVdata(£0 ^ 

^theory 



0-(&)modified = (^theory (^)) V4 if ^data(^) = . (4£ 



Obviously, if theory and experiment match very closely the a's in (|4*Tj) and ()48|) are 
essentially identical. This happens in the bulk of the distribution. For the tails 
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m] 


[1] 


[2] 


[3] 


[5] 


[7] 


[10] 


[20] 


[30] 


N 

1 "eve 


10698 


7106 


10698 


10698 


8828 


10698 


9009 


6243 



Table 1: The number of events N eve for the maturities [m]. 

of the distributions, when iV data is equal to zero or one event, the a (k) modified is 
smaller than o"(fc) natU rai- Hence, more emphasis is put on these tail points by the 
minimization procedure. The tests performed using both forms of cr(fc) confirm this 
argument. The use of (|4*Hj) leads to a very good agreement between theory and data 
in the tails as can be seen in the Figures. 

5 The Federal Reserve System data. Determina- 
tion of iVdata^ ^bin) f° r each allowed maturity 

The FRS [3] data gives in successive working days the daily average interest rate 
/(day) between banks. They are given in Table Qfor the eight following maturities 

{[m]} = {[1], [2], [3], [5], [7], [10], [20], [30]} in years . (49) 

For many of these maturities, the data extends (10698 events) without break from 
1962 to today (we have chosen November 4 as our last event). For a few maturities, 
the data is restricted to one or more sub-periods within the 1962-today period. The 
total number of events for each maturity is given in table (0). In particular, it should 
be stressed that the data for [m] = 20 years consist of two disjoint periods. We have 
not been able to discover if, during theses two periods, the way of computing the 
daily averages are coherent. We have chosen to group the two sets but we will be 
led to caution comments. 

The interest rates are given in percent by a number with exactly two decimal 
figures, i.e. by a integer number in basis points. It allows the definition of the 
following meaningful distributions for a lag L of one day. These distributions have 
to be compared to the corresponding theoretical distributions. 

The empirical N(v, I) and N(v, i) density distributions 

For v and I both expressed in basis points and for each of the allowed maturi- 
ties, the empirical distribution N(v, I) is the number of occurrences when the 
interest rate of the FRS on some day was / and the next day I + v. These 
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numbers are statistically very small. It is useful to consider the compounded 
empirical distribution N(v, i) where v is still expressed in basis points but the 
interest i is expressed in percent by grouping the days when the interest rate 
is almost i 

More precisely, for the data, the N(v, I) (v and / in basis points) density 
distribution is defined as follows 

N(v, I) = number of days when J(day) = / and /(day + 1) = I + v . (50) 

Note that, in order to obtain /(day+1), non working days are simply discarded. 
The average discretised N(v,i) (v in basis points and i in percent) density 
distribution is precisely defined by 

100i+50 

7t(v,i)= Y, N ( V J)- (51) 

J=100i-49 

It is the number of days when the interest rate was between i — 1/2 percent 
( exactly 100* — 49 basis points) and i + 1/2 percent (exactly lOOz + 50 basis 
points) and the interest rate has moved by an amount v basis points by the 
next day. 

The empirical w(I) and w(i) interest distributions 

For each allowed maturity, the empirical interest distribution w(I) is defined 

as 

w{I)= YN(vJ) (52) 

where I is in basis point. It is the number of days when the interest rate was 
/ in basis point. The average discrete version w(i) of (|52~j) is defined by 

100i+50 

w(i)= J2 ^OO- (53) 

r=100i-49 

where i is an integer giving the interest rate in percent, it is simply the number 
of days when the interest rate was between (i — 1/2) and (i + 1/2) percents. 
The total number of events w = X)jW(i) depends on the maturity and is given 
in Table ©. 

The w(i) are given in Table (J2J) for all maturities of [m] years. In preceding 
papers it was realized that bins with less than about one thousand events 
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i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


w(i) {1] 


356 


402 


664 


1269 


1674 


2147 


1110 


1147 


582 


w(i) [2] 


60 


418 


293 


569 


599 


1530 


759 


835 


561 


w(i)[3] 


9 


369 


524 


1246 


1254 


2166 


1642 


1334 


732 







33 


45 


130 


1278 


2010 


1791 


1563 


778 


w(i) [7] 








154 


495 


693 


1705 


1775 


1582 


880 


w(i)[io] 








28 


1485 


1246 


1949 


1874 


1629 


929 


w(i)[2o] 











1020 


1320 


1974 


1389 


1341 


396 


w(i)[ao] 














289 


1076 


968 


1358 


981 



Table 2: The empirical w(i)\ m ] for % — 1 — 9 and all maturities [m]. 

present too much statistical fluctuations. Hence it is useful to group the data 
in bins. We have chosen this bin definition identical for each of the allowed 
maturities. In Table (j3J), the w(ibin) are given, as an example, for maturity 
[m] = 1, together with a definition of the bins. 

The empirical N<i ata (v, ibin) density distributions 

The final empirical density distribution N(v, i^n) which enters the ^-squared 
merit function is taken as 

N data (v, iun) = W W N(v, i) (54) 

6 The definition of iVtheory an d the %-squared merit 
function 

For each maturity, the empirical distribution iVdata(w,i) has to be compared to 
the effective distribution iV thcory (t>, i) computed from the theoretical distribution 
evaluated at lag L = 1 day, which means for I = N L . This distribution p(N L , v, i) is 
computed from the seed by Nl successive convolutions. 
We have 

N t heory{v,i) = W(l) p(N L , V, l) 
^theory 0, «bin) = ^theory ( V, l) ■ (55) 
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i 


10 


11 


12 


13 


14 


15 


16 


17 


18 




469 


176 


202 


105 


160 


131 


65 


38 





w(i) m 


445 


294 


196 


157 


190 


126 


55 


18 







333 


351 


179 


173 


234 


88 


59 


4 







341 


338 


222 


216 


240 


94 


39 








w(i) [7] 


332 


306 


300 


250 


233 


96 


26 










306 


326 


331 


258 


229 


97 


10 








W(l)[20] 


204 


384 


438 


251 


211 


74 


5 








^(*)[30] 


282 


364 


413 


295 


177 


39 












Table 3: The empirical w(i)\ m ] for % = 10 — 18 and all maturities [m]. 



^bin 


1 


2 


3 


4 


5 


6 


7 


8 


% G 2 b in 


1-3 


4 


5 


6 


7 


8 


9-10 


11-17 


W(i hi n)[l] 


1165 


1269 


1674 


2147 


1110 


1147 


1051 


877 



Table 4: Example of the empirical w(iun) for the maturity [to] = [1]. The other [to] 
can be computed from Table (jHJ). 
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The %-squared merit function 

For each maturity [m], the following ^-squared function is defined 

x [m] = L L • (56) 

i bin =l v=-300 v\u,bm) 

In this formula, the numerator is the square of the difference between the 
empirical distribution (J51|) and the theoretical distribution (J55|) for a given 
maturity.. It is a measure of the discrepancy between theory and experiment. 
The denominator is the square of the modified error a(v ,bm ) According to the 
previous discussion around (jUJ), this cr(t>, bm) is taken as 

c(f,ibin)modificd = (iVdata(^,^bin)A'theory(f,ibm)) 1/4 if -Ndatafa, «bin) ^ 
0-(v,ibin)modificd = (theory 0, «bin) ) if iVdatafa, Z b in) = (57) 

With our definitions, the bins where the number of events is zero do not 
contribute to the x-squared and hence can be removed or kept without altering 
the results of the minimisations. 

Finally, let us note that the limits on the summation over v, [—300, 300], have 
been chosen at the point where the distributions in v have become empirically 
completely negligible. For example, increasing the limits from 300 to 400 
basis points does not change the values of the parameters provided by the 
minimization method in any appreciable way. 

7 Seed Parametrisation 

In preceding articles [I], 0, we have shown that the FRS distributions can be fitted 
very closely by using Pade Approximates [0, 4] whose coefficients follow rather simple 
scaling laws (see (jHlJ)). 

7.1 The form of the Seed. Determination of the parameters 



Extrapolating the scaling laws discovered in pQ, [2] for lags equal to an integer 
number of days to values of the lag small compared to one day, as is explained in 
Appendix (jDj), one is led to state the following form for the discretized seed 

1 a(i) 

s L S{vs v ,isi) = — — — d(i] forw^O. (58) 



18 



It depends on three independent parameters a, 7 and d which themselves depend 
on the initial interest rate % and on the maturity [m]. In this equation, using our 
choices of Sl (J35]l . sy (|37j) and sj (I59*|l . we see that v and z are integers, that the 
variation of interest rate V = syv has the dimension of basis point, that the initial 
interest rate I = Sji has the dimension of percent, that the contribution of slS to 
the probability distribution has the dimension of the inverse basis point and that a, 
13 and d have no dimension. 

It is first convenient to define V tra nsiticm as the variation of interest rate V at 
the place where the transition between a constant behavior (in v) and the power 
decrease of the seed occurs. More precisely, it is defined as the value for which 
the second term in the denominator of (J58|) becomes equal to the first term in this 
denominator. Thus when the second term becomes 1. It is linked to the parameter 
7 by 

Vtransition ~ l' 1 • (59) 

Using the definition (J36j) . we see that this corresponds to a value of the interest rate 
variation Vtransitim expressed in basis point of 

^transition ^transition ■ 

(60) 

It is then easy to understand the meaning of the three parameters: 

• The parameter a is linked, in first approximation, to the constant behaviour 
of the seed for (non-zero) variations of the interest rate appreciatively smaller 
than the transition value defined in (|59|). 

• The parameter 7 is linked, in first approximation, to the transition value of 
the variation of the interest rate (J59|) between a constant behaviour (for small 
v) and a power decrease of the seed for large v. 

• The parameter parameter d is linked to the asymptotic power decrease of the 
seed when the variation of the interest rate is large compared to % 'transition- 

Guided by the findings of P, 0, we have taken the following dependence of the 
three parameters in terms of the initiail interest rate: 

a(i) = ai e a ^- io) 

7 (l) = ry ie W(i-i0) 

d(i) = d ie ^- 4 °) + d 3(*-*0) 2 _ ( 61 ) 

Again i is in percent and Iq is arbitrarily but very conveniently chosen to be 6.5 
percents (see the discussion in Appendix (|B|l). At this point, there are thus seven 
parameters which are maturity dependent. 
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7.2 Determination of the parameters 

To calibrate the model, we have then proceeded along the following steps. 

• We have first computed, for each maturity [m] , the values of the seven param- 
eters (q;i, a 2 , Ti) 72, dii d 2 , d 3 ) which minimize the corresponding X[ m y 

• Looking at these minimal values, we have discovered that linear extrapolation 
in [m] are approximately at work. More precisely, the following parametrisa- 
tion with straight lines is suggested 



«i(H) 


= (a n + a 12 ([m 


- Ho)) 10- 


a 2 ([m}) 


= (a 21 + a 22 ([m 


- Ho)) 10" 


7i(H) 


= (7n+7i2(H 


- Ho)) 10" 


72(H) 


= (721 + 722 (H 


- Ho)) 10- 


di ( [m] ) 


= (d u + di 2 ([m] 


- Ho)) 


d 2 ([m]) 


= (d 21 + d 22 ([m] 


- Ho)) 10" 


d 3 ([m}) 


= {dn + d 32 {[m) 


- Ho)) 10- 



The [ m\ which appears in these formulae is arbitrarily chosen to be [15] (years) 
close to the middle of the interval [[m] = [1], [in] = [30]]. The normalisation 
10" 5 , 10" 6 , . . . have been chosen conveniently. Clearly, the parameters with 
a 2 as their second index are related to the slopes of the straight lines while 
the parameters with their first index are the ordinates of the lines for 

[m] = [m]o- Altogether, there are now fourteen constant parameters. 

• To determine the empirical values of the fourteen parameters we have chosen 
to minimize the total merit function 

Xtotal — Xfm] (02) 

[m]e{[m]} 

which is simply the sum of the x- sc L uare d for the eight maturities (|49|). See 
the discussion in Appendix JOJ) where the results are summarized in Tables 

• It turns out that only eight of the fourteen parameters are really relevant and 
that six of them, namely a 2 x, a 22 , di 2 , d 2 i, d 22 , d 3 \ can be safely put to zero (see 
the discussion in Appendix (|D|) 

"21 = "2i = di2 = d 21 = d 22 — d 31 — . (63) 
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a n 


«12 


7n 


712 


721 


722 


d u 


d 3 2 


value 


2.39 


0.018 


5.78 


0.20 


-2.90 


-0.100 


3.030 


0.909 


C.I. 


± 0.08 


± 0.004 


± 0.99 


± 0.04 


± 0.49 


± 0.032 


± 0.024 


± 0.067 



Table 5: The selected a, 7, d parameters with the choice io = 6.5, [m]o = 15 and 
a 2 i = «2i = ^12 = ^21 = ^22 = d 3 i = 0. The C.I. are the confidence intervals. 
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[1] 


[2] 


[3] 


[5] 


[7] 


[10] 


[20] 


[30] 


a x 10 5 


2.14 


2.15 


2.17 


2.21 


2.24 


2.30 


2.4 


2.6 



Table 6: The values of a x 10 function of [ml. 

• A last minimisation, taking again Xtotai as ^ ne m erit function and taking into 
account the restrictions (J6ll|) . has produced the final values of the eight re- 
maining non-zero parameters as given in Table (jSJ). 



7.3 Comparison with the data. Discussion 

We are now in a position for a detailed comparison of the results of our fits with the 
data. The values obtained for the parameters are summarized in Table © for a, in 
Table (JTj) for 7, in Table (jH)) for d and in Table Q for v transition- They are given in 
natural units (see [T] for a precise discussion of these units) . 

The parameter a 

The parameter a does not depend of the initial interest rate and is a slowly 
increasing function of the maturity [m]. Its values (multiplied by 10 5 ) are 
expressed in inverse basis points and given in Table (0). 



The parameter 7 

The parameter 7 depends both on the maturity and on the interest rate. Its 
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7 x 10 6 


i=l 


i=3 


i=5 


i=7 


i=9 


i=ll 


i=13 


[m] = [l] 


6.85 


5.08 


3.77 


2.8 


2.07 


1.54 


1.14 


[m] = [2] 


7.72 


5.61 


4.08 


2.96 


2.15 


1.57 


1.14 


[m] = [3] 


8.66 


6.17 


4.39 


3.13 


2.23 


1.59 


1.13 


[m] = [5] 


10.8 


7.38 


5.05 


3.46 


2.37 


1.62 


1.11 


[m] = [7] 


13.3 


8.74 


5.75 


3.78 


2.49 


1.63 


1.07 


H=[io] 


17.9 


11.1 


6.86 


4.25 


2.63 


1.63 


1.01 


[m]=[20] 


43.9 


22.3 


11.3 


5.71 


2.89 


1.46 


0.742 


[m]=[30] 


98.6 


40.9 


16.9 


7.02 


2.91 


1.2 


0.499 



Table 7: The values of 7 x 10 6 as a function of [m] and i. 

dimension is the inverse of the basis point raised to the power 1/d. Its is given 
in Table (0) multiplied by 10 6 . In fact v tra nsition defined before, directly derived 
from 7 and given in Table 0, carries a more intuitive meaning. 



The parameter d 

The parameter d (which is a pure number) depends on the maturity and on 
the initial rate. It is given in Table (JHJ). It should be remarked that it is 
empirically always close to 3. One may even wonder if the exact value of 3 
would not be a critical exponent of the problem and have general validity. 



The transition point v transition 

With the values of the parameters in (j5J, the values of the transition points 
Vtransition are given in Table © for the odd z's and all the maturities. One sees 
that, as expected by intuitive arguments, v trans i tion is an increasing function 
of the initial interest rate and a decreasing function of the maturity. 



The x- sc L uare d for the best parameters 

The total x-squared as well as the x-squared for all maturities can be found in 
table (JT2J). In general, as discussed in the appendix, the reduced %-squared are 
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d 


i=l 


i=3 


i=5 


i=7 


i=9 


i=ll 


i=13 


H=[i] 


2.92 


2.98 


3.02 


3.03 


3.01 


2.95 


2.87 




m] = [2] 


2.92 


2.99 


3.02 


3.03 


3.01 


2.96 


2.88 




m] = [3] 


2.93 


2.99 


3.02 


3.03 


3.01 


2.96 


2.89 




m] = [5] 


2.95 


3.0 


3.02 


3.03 


3.01 


2.98 


2.92 




2.96 


3.0 


3.03 


3.03 


3.02 


2.99 


2.94 


[m]=[10] 


2.99 


3.01 


3.03 


3.03 


3.02 


3.0 


2.97 


[m]=[20] 


3.07 


3.05 


3.03 


3.03 


3.04 


3.06 


3.09 


[m]=[30] 


3.16 


3.08 


3.04 


3.03 


3.06 


3.12 


3.21 



Table 8: The values of d SIS db function of [m] and i. 



^transition 


i=l 


i=3 


i=5 


i=7 


i=9 


i=ll 


i=13 




m 


= [1] 


59 


60 


62 


68 


78 


93 


117 




m 


= [2] 


56 


57 


61 


67 


77 


92 


115 


[m]=[3] 


53 


55 


59 


66 


76 


90 


113 


[m] = [5] 


48 


52 


56 


64 


74 


88 


110 




M=P1 


44 


48 


54 


62 


72 


87 


107 


[m]=[10] 


39 


44 


51 


59 


70 


85 


104 




m] = [20] 


26 


34 


43 


54 


67 


81 


97 


[m]=[30] 


19 


27 


37 


50 


65 


80 


92 



Table 9: The values of v tr ansition (see (jH^j) . (??)). as a function of [m] and i. 
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close to one and hence very good, except for [m] = 20 where there are some 
uncertainties in the data and to a lesser extend for [m] = 5. This is reflected 
in the quality of the agreement between the empirical distributions and the 
data as seen in the plots below. 

7.4 A few plots 

All the plots, for lag one and higher, for all bins and for all maturities show good to 
very good agreement between our theoretical curves and the data provided by the 
FRS site jsJj. Both the central region where the variation v of the interest rates is 
around zero and and the tails where v becomes large are very well accounted for. We 
have selected randomly eleven plots (fT|)- (jll|) as examples of the fits. In particular, 
the domains of the selected plots are given in detail in the figure captions. 

8 Conclusions 

In this paper, a microscopic theory of the term structure of interest rates has been 
developed. Convolution techniques, implying about ten successive convolutions, 
combined with time translation invariance lead to a time scaled theory where the 
term structure for practical lags (one day or more) can be deduced from a seed 
function living at a relatively small lag scale (a few seconds). 

The previously discovered scaling laws, which were found to be valid at lags 
expressed in days, suggest forms for the possible seeds and imply a discretisation of 
the problem. 

Our new results show that the scaling law assumptions are even simpler at the 
microscopic lag scale. Indeed, it is shown that the FRS data are amazingly well 
reproduced (for a lag of one day or two days but the results easily extend to higher 
lags) by assuming that the seed has the critical form of a self-organized [5| econo- 
physical system . In other words, a very simple power law behavior emerges with 
essentially only one v~ d term. The exponent d is of the order of three in close 
agreement with the tail behaviours obtained previously using the Hill ^Uj estimator. 

These results open the door to two major issues, one rather theoretical and the 
other more practical and pragmatic. 

Since the seed has such a simple scaling form, it suggests the existence of an 
underlying statistical model. The discovery of the basic ingredients and laws leading 
in a natural and systematic way to the scaling would be a major achievement which, 
if attained, may also lead to new theoretical insights in related but different contexts. 
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We expect that a self-organized structure may be at work, with d — 3, or close to 
three, as critical exponent for all the maturities. We hope to come back to this issue. 

Finally, one may wonder how this theory and its predictions can be efficiently 
used in the context of risk management (e.g. Value-at-Risk computation) and fixed- 
income derivative pricing. For the time being, studies are under way to measure the 
add-value of using this model instead of the traditional ones. First results will be 
provided in a near future. 



25 



A A naive introduction to Dirac type distribu- 
tions 



A.l Formal definition 

A distribution G(V) is a continuous linear functional on some space of functions 
of a real variable V. The Dirac S(V) distribution, which has compact support, is 
defined on a continuous C°° function /(V) at V = by 

(5(V),f(V)) = f(0) . (64) 

The jth derivative of a distribution of compact support G(V) is defined as 

(G®{V), f(V)) = (-iy (G(V), fW(V)) (65) 



where 

M(v\ = - 



f^V) = ^g 1 (66) 



is the jth derivative of the function /(V). 

In particular, the second derivative of the Dirac distribution 

(5®(V)J(V)) = f®(0) (67) 

is the second derivative of /(V) evaluated at V = 0. 

Without being mathematically rigorous, these definitions can be seen in a more 
intuitive way. As Laurent Schwartz, the inventor of the general distribution concept, 
used to caution : The formulae given below, with physicist notations, have to be 
used with great care and have to be justified in a detailed and precise way. But he 
also recognized that they were very convenient to guess properties of distributions. 



A. 2 Discussion and approximations for the Dirac distribu- 
tions 

The distribution S(V) is a generalized function which can be thought as being zero 
for V < and V > and infinite at V = in such a way that, in formal agreement 
with (JUIJ) 

/+oo 
S(V) f(V) dV = /(0) . (68) 
-oo 

There are many functions which can approximate the delta function. Let us give 
two. 
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1. The step-type function 

rD(V,s) = for V< 



8(V) = km D(V,s) < 

s— >0 



D(V, s) — - for -~<F<§ (69) 



lD(V,s) = for §<v 

has compact support. It is not difficult do the Riemann integration and prove 
that for some V 

-\<V< S - (70) 

one has 

J + ™ D(V,s)f(V)dV = f(V)(^J + ~ 2 D(V,s)dV ) j 

= f(V). (71) 
The limit for s — > reproduces clearly (pj%j) as f(V) — > f(0). 
2. The Dirac distribution can also be approximated by the continuous function 

5(V) = hm£(V,A) 

E ^ A) = ma^T (72) 

This form is very interesting, as it shows that the S(V) function is the boundary 
value of a continuous function of two variables. For positive A, this function 
is purely positive for all V. 

A step-type approximation for the second derivative of the delta function is 
given as follows. Consider the intervals A = [—3s/2,—s/2],B = [— s/2, +s/2], A = 
l+s/2, +3s/2], the step function which is zero outside the three intervals, has value 
1/s in the intervals A and C and — 2/s in the interval B and let s —>■ 0. 

B Formal discussion of the Master Equation. Con- 
tinuous variables 

In the following it is shown that the formal solution of the master equation is fully 
compatible with the convolution equation ()18|) whatever be the intermediate time 



27 



Suppose that, the probability distribution has a Taylor power expansion, around 
L = 0, of the form 

oo j n 

p(L,y,/o) = E -yP^iVJo) . (73) 



n=0 n - 



where p( n \V,Io) is the n'th order derivative of p with respect to L evaluated at 
L = 

p^(V,I ) = d n p(L,V,I ) . (74) 

L L=0 

The zeroth order p'°'(V, Jo) is known by Eg. (|21|) and the firth order pW(V, Jo) by 
Eq.pijl 

p(°H^/o) = p(0,V,I o ) = <W 

p (1) (V,J ) = S(V,I ). (75) 

Introducing this power expansion in the master equation (|27j) and equating, in the 
left and right hand sides, the coefficients of the same powers in L, one gets the 
recurrence equations 

/+oo 
pfr-Vfa-IiJjSili-IoJddli, n=l,...,oc. (76) 
-oo 

It is easily seen that the equation for n = 1 is an identity. For n = 2, one finds 

p®(I f - J , h) = / S(I f - h, h) S(h - J , I ) dh (77) 

J — oo 

where the dummy integration variable has been called I\. For arbitrary n > 1, one 
obtains 

r+oo r+oo r+oo 

p(»+V(I f -I Q ,I ) = / dlj dl n ^... dh 

J —oo J —oo J —oo 

S(I f - I n , I n ) S{I n - I n _ x ) . . . - Jo, /o) 

/•-foo ( n \ ( n \ 

= \U di pj (n^+i-^^)! > .(78) 

It is not difficult to check that the formal power expansion (JZHJ), f7H|) . (J77|) . (J7H|) 

satisfies the convolution equation (fTH|) without any further condition. Hence, it can 
reasonably be supposed that one can focus on the master equation (|2*7|) together 
with its natural boundary conditions. 
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C Master Equation with a simplified Scaling Law 
and no Iq dependence. The Gauss distribution 
solution 

As suggested in section (J5J), it is convenient to study analytically the solution of the 
master equation with the seed (|2*9"j) 

S{V,h) = K{I Q )d 2 v 8{V) . (79) 

In order to simplify the problem let us limit ourselves to suppose that the dis- 
tribution p does not depend on I 

p(L,V,I )=p(L,V) (80) 

as well as the seed ()29|) and hence 

k(/o) = « • (81) 

Eq. pUj) becomes 

d L p(L,V) = Zd 2 v p(L,V) . (82) 

To simplify the problem even further, "scaling laws" which are approximately 
but not exactly verified by the data j2], are assumed. It has to be emphasized that 
the data do not support exactly these scaling laws. Hence the results we will obtain 
in this section cannot be hoped to be correct. On the other hand, it is worth to 
quote them since they have as a consequence the very common but wrong belief 
that, in a natural way, distributions should fall off at large V as exponentials. 

The ultra simplification of the empirical results found in pQ, j2] is equivalent 
to the statement that the distributions depend on the reduced variable V (J83|) an d 
scales (J83J) in the natural way as 1/ \fh 

p(L,V,I ) = 5^ (83) 

9 = 7z- (84) 

Replacing p by its guess (j8Hj) in (|3Uj) . one finds the final equation 

2Kd 2 ~p{V) + (vp(V)^ = . (85) 



29 



The solutions of this equation can easily be found. Indeed, it can be written 
whose general solution is 



dy l2ndy p(V) + V p(V) ) = (86) 



2Kdyp{V)+Vp{V) = R 1 (87) 

where Ri is an arbitrary constant. To solve this first order differential equations, 
the homogeneous equation is solved 

2Kdyp h (V) + Vp h (V) = (88) 
giving the homogeneous solution p h (V) 

p h (V) = He-£. (89) 

Assuming (the method of variation of constants) that H depends on V, introduce 
in (J%7|) to obtain the equation 

d v H(V) = ^ (90) 
which can be solved by a simple integration 

H(V) = ^ [ V dx + R 2 (91) 



2k Jo 

where R% is the second constant of integration. One has now to plot the final solution 
obtained by replacing the solution for H fl9*T|) in the homogeneous solution (|%9")l 

p(V) = H(V)e~£ 

/ dxe^ +R 2 )e~ . (92) 



.2k Jo 

In fact, we are looking for a distribution function p(L, V, Jo) which is, in first 
order symmetric in V, hence for a scaled function p(V) which is also symmetrical 
in V. Since by (J89|) the homogeneous solution p (V) is symmetrical, one looks for a 
symmetrical H(V) solution. This clearly implies that Ri should be zero. 

We thus conclude that the solution of the master equation with strict scaling law 
leads to the familiar Gauss type behavior 

p(V) = R 2 e-£ (93) 

which is not at all sustained by the data as the main result of our preceding inves- 
tigation has shown that the asymptotically the term structure decreases as a power 
of 1/V d with d around three to four. 
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D Form of the Seed. Phenomenological Discus- 
sion 



In preceding articles PP, 0, distributions fitting very closely the FRS data have been 
obtained using Pade Approximants [0,4] (see ()94|l) i.e. a polynomial of zero degree 
in v in the numerator divided by a polynomial of fourth degree in the denominator 

a(i) , . 

P(v,t)= Y- 4 ■ 94 

Moreover, it was shown that the Pade coefficients (a(i), b(i), c(i)), which also depend 
on the maturity, follow rather simple scaling laws. Extrapolating these scaling laws 
for values of the lag small compared to one day have lead us to guess suitable forms 
for the seed. Though, the asymptotic behavior | v \ ~ d for large v with d equal four 
which follows from (|94|) is not incompatible with the data, Hill estimators [IQ] were 
pointing towards a somewhat smaller value of d. Precise values can be found in pQ. 

Guided by this work and by general ideas about scaling laws and self-organized 
criticality, we have progressively been led to a simpler educated guesses for the seed, 
namely ([58)1 . where three parameters a, 7 and d, in principle, depend on the initial 
interest rate % and maturity [m]. For the dependence in i, we were led to ()6ip where 
Zo can be arbitrarily chosen. The seven coefficients themselves are thought to be, in 
first approximation linear in the maturity m. More precisely, we have thus used 
and (|62|) with their convenient normalisations. They depend on initially on fourteen 
independent constant parameters. 

At first, we have chosen, as is generally advisable to diminish the relationships 
between the relevant parameters, the arbitrary iq an [m]o in the middle of their 
respective domains and more precisely io = 6 and [m] = [15]. Starting with this 
form of the seed and allowing for Nz = 10 iterations (see the discussion following 
eq.(J3U), after minimizing the total x-squared we find the results of Tables (fTUj) and 
(fTTj) for the fourteen parameters in (JB2j). The values of the resulting xf m ] an d Xtotai 
are given in Table ()12|) together with the number of degrees of freedom Nf and the 
reduced x- sc i uare d X 2 /Nf for each maturity. 

Let us now discuss the results of the Tables (fT2j). (fTU|) . (fTTj) and their implications 
for the choice of the parameters. 

The x 2 

From Table (fT2|) we see that the reduced x- sc L uare d are rather good except 
for [m] = [20]. Let us recall that there is some incoherences for this maturity 
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"11 


"12 


"21 


"22 


7n 


712 


721 


722 


value 


2.46 


0.022 


0.005 


-0.015 


6.53 


0.230 


-2.90 


-0.119 


C.I. 


± 0.08 


± 0.0042 


± 1.57 


± 0.122 


± 0.99 


± 0.039 


± 0.49 


± 0.032 



Table 10: The a and 7 parameters for the initial arbitrary choice of i$ = 6 and 
[m] = [15]. The C.I. are the confidence intervals. 





dii 


di2 


^21 


d 2 2 


d 3 i 


d 3 2 


value 


3.015 


0.0001 


0.04 


0.005 


-0.01 


0.907 


C.I. 


± 0.024 


± 0.0005 


± 2.47 


± 0.161 


± 0.82 


± 0.067 



Table 11: The d parameters for the initial arbitrary choice of io = 6 and [m]o = [15] 
. The C.I. are the confidence intervals. 





[1] 


[2] 


[3] 


[5] 


[7] 


[10] 


[20] 


[30] 


total 


x 2 


1064 


1164 


930 


1440 


1090 


996 


1832 


828 


9347.44 


N f 


1099 


931 


959 


871 


853 


793 


691 


651 


6911 


X 2 /N } 


0.97 


1.25 


0.97 


1.65 


1.28 


1.26 


2.65 


1.27 


1.35 



Table 12: The x 2 , the reduced x 2 '/Nf and the number of degree of freedom Nf 
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as the data is composed of two disjointed sets of points corresponding to two 
disjoined periods in time. 

The parameters a 

From Table (jl(J|) . we see that two parameters «2i and a 22 are very close to 
zero with errors much larger than their values. This obviously points to a zero 
value for the combined parameter a 2 which, we recall, sets the dependence in 
i of the parameter a. It means that whatever be the level of the interest rate, 
for very short intervals the probability of the occurrence of a certain small 
variation (before the decrease in \v \~ d plays a role) does not depend on the 
interest rate. This is an interesting conclusion. From Table (jlUj) . we also learn 
that the two other parameters (an and 0:12) are meaningful. We conclude 
that the parameter a, which essentially does not depend on the initial interest 
rate, increases appreciably with the maturity. 

The parameters 7 

The four 7 parameters are all significant. This means that % 'transition ((EH) 
where the constant behaviour of the seed tends to become a power behaviour is 
highly dependent on the maturity and on the initial interest rate. As discussed 
in the main part of the text, this behaviour conforms to our expectations. 

The parameters d 

From Table (jTTjl. we see that of the six d parameters, four, namely d± 2 , d 22 , d 22 , ^31, 
can safely be put to zero. This means that d is almost constant for every ma- 
turity and every initial interest rate. Except for du which as expected is very 
close to 3, the only parameter which survive is d% 2 . It follows that d has a 
joint behaviour linear in maturity and parabolic in the initial interest rate. 
Amazingly enough the arbitrary parameters io = 6 and [m] = [15] which 
were chosen in the middle of their respective domains were almost optimal 
choices which are further discussed below. 

In order to obtain the best value of zo leading to the best adjusted parabola 
for the dependence of d in i, we have rerun our minimization programs putting to 
zero the six un-significant parameters a 2 i, a 22 , di 2 , d 2 %, d 22 , d 3 i but allowing the best 
choice for the arbitrary parameter i . The results are given in Table (JHJ). They show 
that i around 6.5 would have even be a better choice to make di 2 , d 2 i, d 22 , d%\ equal 
to zero. In the main part of the text we have adhered to this choice. 

Acknowledgment 



33 



The work of Jean Nuyts was supported by the Belgian F.N.R.S. (Fonds National 
de la Recherche Scientifique) and the work of Thomas Alderweireld partly by the 
Belgian I. I.S.N (Institut Inter-universitaire des Sciences Nucleaires) and the DEXIA 
group. 



34 



References 

[1] Nuyts, J. and Platten, L, Physica A, 299, 528-546 (2001). 

[2] Alderweireld, T. and Nuyts, J., Physica A, 331, 602-616 (2003). 

[3] Federal Reserve Statistics. Historical data. 

http://www.federalreserve.gov/releases/hl5/data.htm. 

[4] Bharucha-Reid, A.T., Elements of the theory of Markov process and their Appli- 
cations, New- York, McGraw-Hill (1960). 

Papoulis, A., Probability, Random Variables and Stochastic Processes, New- York, 
McGraw-Hill (1984). 

[5] Levy, P., Theorie de I Addition des Variables Aleatoires, Gauthier-Villars (1937). 

[6] Mandelbrot, B.B., Journal of Business, 36, 394-419 (1963). 

[7] Mandelbrot, B.B., Fractals and Scaling in Finance, Springer (1997). 

[8] Weron, R., Int. J. Modern Physics C, 12, 209-223 (2001). 

[9] Back, P., Tang, C, Wiesenfeld, K., Phys. Rev. Let, 59, 381-384 (1987). 

[10] Hill, B., Annals of Statistics, 35, 1163-1173 (1975). 

Nuyts, J., Inference about the tail of a distribution. Improvement of the Hill 
estimator, to be published 



35 



Figure Caption 



Figure 1. The comparison between the theoretical curve (line) and the empirical 
curve (crosses) for [m] = [1] and for the first bin %u n = 1 which contains the 
initial interest rates i = 1,2,3. The related total number of events is found 
in (0J). The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors 
equal approximatively to the square root of the number of events. The curves 
for the bins iu n = 4 and iu n = 8 and for the same maturity [in] = [1] are given 
in Figures (0) and (j3J) respectively. The analogous curves for the other bins 
present fits of equivalent quality. 

Figure 2. 

The comparison between the theoretical curve (line) and the empirical curve 
(crosses) for [m] = [1] and for the bin i^n — 4 which contains the initial interest 
rate i = 6. The related total number of events is found in (@J). The horizontal 
axis is in basis points. The vertical axis represents the number of events. The 
size of the arms of the crosses are estimated errors equal approximatively to 
the square root of the number of events. The curves for the bins %u n = 1 and 
hin = 8 an d for the same maturity [m] = [1] are given in Figures and Q 
respectively. The analogous curves for the other bins present fits of equivalent 
quality. 

Figure 3. 

The comparison between the theoretical curve (line) and the empirical curve 
(crosses) for [m] = [1] and for the bin iu n = 8 which contains the initial 
interest rates i = 11 — 17. The related total number of events is found in 
P)). The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors 
equal approximatively to the square root of the number of events. The curves 
for the bins iun = 1 and iun = 4 and for the same maturity [m] = [1] are given 
in Figures (JTJ) and (0) respectively. The analogous curves for the other bins 
present fits of equivalent quality. 

Figure 4. 

The comparison between the theoretical curve (line) and the empirical curve 
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(crosses) for [m] = [3] and for the bin = 3 which contains the initial 
interest rate i = 5. The related total number of events can be read in Table 
P)). The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors equal 
approximatively to the square root of the number of events. The analogous 
curves for the same maturity [in] = [3] for the other bins, as well as those for 
maturity [m] = [2], present fits of equivalent quality. 

Figure 5. 

The comparison between the theoretical curve (line) compared with the em- 
pirical curve (crosses) for [m] = [5] and for the first bin iu n = 4 which contains 
the initial interest rate i = 6. The related total number of events can be read 
in ©. The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors equal 
approximatively to the square root of the number of events. The analogous 
curves for the same maturity [in] = [5] for the other bins, as well as those for 
maturity [m] = [7], present fits of equivalent quality. 

Figure 6. 

The comparison between the theoretical curve (line) compared with the empir- 
ical curve (crosses) for [m] = [10] and for the first bin iu n = 5 which contains 
the initial interest rate i = 7. The related total number of events is found in 
(j2*j). The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors 
equal approximatively to the square root of the number of events. The analo- 
gous curves for the same maturity [in] = [10] for the other bins present fits of 
equivalent quality. 

Figure 7. 

The comparison between the theoretical curve (line) compared with the empir- 
ical curve (crosses) for [m] = [20] and for the first bin % in = 3 which contains 
the initial interest rate i = 5. The related total number of events is found in 
(j2*j). The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors 
equal approximatively to the square root of the number of events. The analo- 
gous curves for the same maturity [in] = [20] for the other bins present fits of 
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equivalent quality. 



Figure 8. 

The comparison between the theoretical curve (line) compared with the empir- 
ical curve (crosses) for [m] = [30] and for the first bin 4m = 7 which contains 
the initial interest rates % — 9— 10. The related total number of events is found 
in @, ©. The horizontal axis is in basis points. The vertical axis represents 
the number of events. The size of the arms of the crosses are estimated errors 
equal approximatively to the square root of the number of events. The analo- 
gous curves for the same maturity [m] = [30] for the other bins present fits of 
equivalent quality. 

Figure 9. 

The comparison between the theoretical curve (line) compared with the em- 
pirical curve (crosses) for the tail (i.e. for v greater than 5 basis points) for 
[m] = [1] for the bin if, in = 7 which contains the initial interest rate i = 9 — 10. 
The horizontal axis is in basis points. The vertical axis represents the num- 
ber of events. The size of the arms of the crosses are estimated errors equal 
approximatively to the square root of the number of events. The analogous 
curves for the other bins and the other maturities present fits of equivalent 
quality. The analogous curves for essentially all maturities all bins, when they 
make sense, present fits of equivalent quality. 

Figure 10. 

The comparison between the theoretical curve for [m] = [5] and for a Lag 
of 2 days (line) compared with the empirical curve (crosses) for the first bin 
hin = 3 which contains the initial interest rates i = 5. The related total 
number of events is found in ©. The horizontal axis is in basis points. The 
vertical axis represents the number of events. The size of the arms of the 
crosses are estimated errors equal approximatively to the square root of the 
number of events. We have checked that the analogous curves obtained by 
further convolutions, for all maturities present fits of equivalent quality. 

Figure 11. 

The comparison between the theoretical curve for ]m] = [1] for a Lag of 2 days 
(line) compared with the empirical curve (crosses) for the tail (i.e. for v greater 
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than 5 basis points) for the bin i bin = 5 which contains the initial interest rates 
% = 7. The horizontal axis is in basis points. The vertical axis represents the 
number of events. The size of the arms of the crosses are estimated errors equal 
approximatively to the square root of the number of events. The analogous 
curves of the other bins present fits of equivalent quality. 
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Figure 1 

Maturity=[l], Initial Interest Rate ~ 1—3, Lag= 
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Figure 2 

Maturity=[l], Initial Interest Rate ~ 6, Lag= 
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Figure 3 

Maturity=[l], Initial Interest Rates ~ 11—17, Lag 
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Figure 4 

Maturity=[3], Initial Interest Rate ~ 5, Lag=l 
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Figure 5 

Maturity=[5], Initial Interest Rate ~ 6, Lag=l 
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Figure 6 

Maturity=[10], Initial Interest Rate ~ 7, Lag=l 
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Figure 7 

Maturity=[20], Initial Interest Rate ~ 5, Lag= 
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Figure 8 

Maturity= [30] , Initial Interest Rates ~ 9—10, La 
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Figure 9 

Maturity=[l], Initial Interest Rates ~ 9—10, Lag=l, Tail 
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Figure 10 

Maturity=[5], Initial Interest Rate ~ 5, Lag=2 
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Figure 11 

Maturity=[l], Initial Interest Rate ~ 7, Lag=2, Tail 
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